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ABSTRACT 

The vaporization time of large drops of ordinary fluids such as water, 
ethanol, or benzene have been correlated by conventional laminar film 
boiling theory. However, liquid nitrogen data has been found to have a 
30 percent shorter vaporization time than predicted by theory. Experi- 
mental and theoretical attempts are made herein to resolve this discrep- 
ancy. 

An asymptotic expression for the vaporization time of large drops in 
film boiling has been derived. This expression seems to correlate the 
nitrogen data. Also, a new expression for the modified latent heat of va- 
porization is derived. This expression is useful for large temperature 
gradients occurring within the vapor film supporting the drop or for fluids 
having small latent heats such as would occur for cryogens or ordinary 
fluids near the thermodynamic critical point. 

NOMENCLATURE 

A area 

* 2 
A dimensionless area, see table I, A/L 

Cp specific heat at constant pressure 

f radiation correction factor (given in ref. 5) 

g acceleration of gravity 

*NASA Lewis. 

fOld Dominion University. 

^Cleveland State University. 
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V 
X 


gravitational constant 
heat transfer coefficient 

dimensionless heat transfer coefficient (see table I) 
heat transfer coefficient from flat plate 
dimensionless % 
thermal conductivity 

[°8c /(p l ‘ P v >e] 1/2 

drop thickness (see table 1) 

dimensionless drop thickness (see table I), 1 / L 

radial coordinate 

temperature 



plate temperature 
saturation temperature 
total vaporization time 

dimensionless vaporization time, see table I 
drop volume 

3 

dimensionless drop volume, see table I, V/L 
axial coordinate 
thermal diffusivity 

dimensionless gap thickness, eq. (15) 
gap thickness 

dimensionless gap thickness, see table I 
dummy variable 
latent heat of vaporization 
modified latent heat of vaporization 
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modified latent heat of vaporization for calculation of vapor gap 
thickness 


( 1 ) A + = A (1 + — — 2. 

20 X 


C AT' 


from reference 9 for small values 


of (C AT)/X 

C AT 
(2) X + = -£ 


-1/3 


= X ( — for large C AT A from 


present paper 

approximate expression for A* for (C AT A) ^ 2 

P 

least upper bound of X* 

approximation to X* for small C 1ST A 

viscosity 

liquid density 

vapor density 

surface tension 

axial velocity 


INTRODUCTION 


Film boiling of liquid drops, called Leidenfrost boiling, plays an im- 
portant role in energy transfer in cryogenic systems. Leidenfrost boiling 
occurs when a liquid drop rests on a cushion of its own vapor resulting 
from heat transfer from a hot supporting surface. For a listing of recent 
publication and more details on the physics of Leidenfrost boiling, the 
reader may refer to references 1, 2, or 3 which contain a comprehensive 
literature summary of Leidenfrost boiling. 

A recent paper by Forslund and Rohsenow (ref. 4) recognized the im- 
portance of Leidenfrost boiling in either nuclear or conventionally fired 
boiler tubes. They used the correlations developed by Baumeister, Hamill, 
and Schoessow (ref. 5) for predicting the heat transfer coefficients in the 
Leidenfrost boiling portion of the boiler tube. The correlations of refer- 
ence 5 for the dimensionless average drop thickness l*, area A*, heat 
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transfer coefficient to the drop h*, and total vaporization t* are shown 
in table I of this report and will be referred to periodically. All these 
parameters are a function of the initial dimensionless volume of liquid 
V* placed on the hot surface. 

The dimensionless vaporization time t* in table I represents the 
dimensionless time for a drop of dimensionless volume V to com- 
pletely vaporize. The equations in table I have been found to correlate 
conventional fluids such as water, ethanol, benzene, or carbon tetra- 
chloride. However, Keshock and Bell (ref. 1) showed that liquid nitro- 
gen Leidenfrost data fell 30 percent below the theoretical curve, although 
the data did follow the general trend predicted by the theory. 

Keshock and Bell (ref. 1) showed that the actual heat transfer coeffi- 
cient was about 30 percent higher than that predicted by the equations in 
table I. This accounts for the shorter vaporization times. They attrib- 
uted this increase in heat transfer rate to vapor instabilities. At large 
dimensionless liquid volumes (V* > 155), the vapor beneath the liquid 
drop breaks through the surface of the liquid. This effect is not con- 
sidered in the theory and could account for the shift downward of the data. 

Baumeister (ref. 6) suggested a second possible explanation for ni- 
trogen data falling below the theoretical curve. In the theory used to 
develop the equations in table I, an asymptotic expansion is introduced 
to account for the convection terms in the energy equation. This expan- 
sion is good only under the condition that 


C (T - 
P P 


T sat> 


« 1 


( 1 ) 


The parameter X is the latent heat of vaporization and (T - T ga p is the 
plate temperature minus the saturation temperature, while Cp is the 
specific heat at constant pressure. The expansion leads to a formula for 
a modified latent heat of vaporization of the form 

i-3 

( 2 ) 


^small * 




C P (T P ~ T saP 


20 X 

The asymptotic expansion used to develop equation (2) is applicable to 
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ordinary fluids since X is generally greater than the product of 
c p (t p - Tgaf)- However, nitrogen has a much smaller value of the latent 
heat than other fluids. Consequently, the expression for the modified la- 
tent heat of vaporization could be in error since equation (1) is no longer 
satisfied, and a corrected value might improve the agreement between ex- 
periment and theory. 

Furthermore, at pressures approaching the critical pressure, equa- 
tion (2) may lead to large errors since the ratio of Cp(Tp - T ga j.) A ap- 
proaches infinity. The energy equation should be investigated to see what 
effect large values of Cp(Tp “ T sa pA have on the expression for X*. 

The purpose of this paper is as follows: 

(1) First, to experimentally investigate the hypothesis put forward by 
Keshock and Bell concerning the effect of vapor instabilities on the total 
vaporization time of large sheets of liquid undergoing film boiling. Va- 
porization time data on both cryogenic and conventional fluids with ex- 
tremely large volumes will be taken and compared. 

(2) Secondly, a least upper bound for X* will be derived as well as 
an exact numerical solution for X *. 


VAPOR INSTABILITIES 


The total vaporization time for a drop is that time required for the 
entire volume of liquid which is placed on a heating surface to vaporize 
completely. The vaporization times of large water, ethanol, and nitrogen 
drops were measured on a 20. 32 cm diameter -1.9 cm thick brass plate 
which was instrumented with chromel-alumel thermocouples. The maxi- 
mum size of drop investigated was 100 ml, which appreciably extends the 
range of the largest volumes investigated previously (i. e. , 30 ml by 
Borishansky (ref. 7)). The plate was heated by a rheostat controlled hot 
plate. 

A plot of the experimental data is shown in figure 1. In this figure, 
the dimensionless vaporization time t* is plotted against the dimension- 
less initial volume V*. The triangular data points are for nitrogen taken 
in reference 1 plus a few additional points (marked with a diamond) taken 
in the present study. As can be seen in figure 1, the dimensionless volume 
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runs from 100 to over 34 000, appreciably extending the range of dimen- 
sionless volumes investigated in reference 1. The actual volume of the 

o o 

nitrogen runs from 0. 16 cm to 40 cm , The extremely large dimension- 
less volumes for nitrogen occurs because of its small value of surface 
tension in the property group that makes up the dimensionless volume V*. 
From the experiments and theory discussed in reference 8, vapor break- 
through can occur for V* greater than 100, that is, the entire range of 
dimensionless volumes shown in figure 1. However, breakthrough will 
be more pronounced for dimensionless volumes greater than 1000. The 
dimensionless volume V* is a logical parameter to correlate vapor 
breakthrough since the property group used to dimens ionalize the actual 
volume is proportional to the cube of the critical wavelength from hydro- 
dynamic stability theory. In reference 8, breakthrough has been shown 
to be related to the critical wavelength. 

In the present experimental study, the vaporization times of water 
and ethanol were measured with dimensionless volumes comparable to 
that of nitrogen. However, because of the relatively high value of sur- 
face tension of water and ethanol (as compared to nitrogen), actual liquid 
volumes from 25 to 100 cc were used. The basic purpose of the present 
experiment was to see if the data on the ordinary fluids would group to- 
gether with the nitrogen data on a V* - t* plot. Since the departure of 
the nitrogen data from the theoretical predictions of reference 5 was hy- 
pothesized in reference 1 to occur because of vapor breakthrough effects, 
a grouping of data would tend to lend further credibility to such a hypoth- 
esis. 

As can be seen in figure 1, however, the ordinary fluids do not group 
with the nitrogen data. The increased amount of vapor breakthrough oc- 
curring in extended masses of ordinary fluids does not appear to appre- 
ciably increase the rate of heat transfer to those drops. Consequently, 
vapor breakthrough may not be the reason for the shorter vaporization 
times of the nitrogen drops. 

LIMITING CASE FOR INFINITE DROPS 

In this section, we will make a theoretical estimate of the effect of 
vapor breakthrough on film boiling of extremely large drops. First, we 
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notice a peculiarity of the expression for the heat transfer coefficient of 
the extended drop model. From table I, 

h* = -L P 4 v* > 155 (3) 

y*l/4 


The derivation of this equation given in references 5 and 9 assumes a flat 
disk liquid drop with a uniform vapor gap beneath it, as pictured in the 
figure in the upper left corner of figure 1. However, for large volumes 
the drop will look more like the lower right hand schematic in figure 1, 
similar to conventional pool film boiling, except that the fluid dept is not 
appreciable. 

As can be seen from equation (3), as V* approaches infinity, the 
expression for the heat transfer coefficient will approach zero. Intui- 
tively, however, we should expect the heat transfer coefficient to be re- 
lated to conventional pool film boiling from a flat surface. The expres- 
sion for pool film boiling heat transfer coefficient h^ is given in refer- 
ence 10 as follows: 



0.410 


f7 3 * ~i^/ 4 

r k x*(p L - p v )p v g 


AT/iL 


( 4 ) 


Nondimens ionalizing equation (4) using the property group in table I gives 
(neglect difference between \* in equation (4) and that used in equation (2)) 


hjp =0.410 

Taking the ratio of equations (3) and (5) gives 


( 5 ) 


^2= 0.25 V* 1 / 4 
h* 


(6) 


Consequently, the break-even point (h^ = h*) occurs at V* equal to 256. 
For dimensionless volumes greater than 256, equation (4) will give a 
higher value for the heat transfer coefficient. It must be pointed out, that 
equation (4) will be only an approximation to the actual value of the heat 
transfer coefficient, since pool film boiling (to which eq. (4) applies) im- 
plies a substantial liquid depth, whereas large unconstrained liquid masses 
have an asymptotic thickness (depth) on the order of 0. 2 inch at most. 
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Nevertheless, let us use equation (5) in place of equation (3) in the 
derivation for the dimensionless vaporization time t when V is 
greater than 256. Following the detailed procedure given in reference 5, 
the vaporization time equation can be expressed as follows: 

t* = 4,5 lnf-^\ for V* > 256 (7) 

\14. 2/ 

This equation is represented by the dotted line in figure 1. 

Intuitively, as the size of the liquid drop increases, the boiling model 
in the upper left hand corner of figure 1 will be transformed to the lower 
right hand model in figure 1. Consequently, one might expect the data 
points to shift from the upper to the lower curve. Perhaps the point where 
the shift occurs is different for different fluids. If this were true, the ex- 
planation given by Keshock and Bell (ref. 1), relative to nitrogen data at 
least, may be valid. 

From high speed movie studies, the vapor domes in nitrogen apppear 
to break open and spew out the enclosed vapor. However, in the water 
drops the domes sometimes remain intact much like a balloon. Perhaps, 
the higher surface tension of water relative to nitrogen accounts for the 
spread between the nitrogen and water data in figure 1. Thus, surface 
tension effects could cause a shift in the boiling model, as will be discussed 
in the discussion section of this report. 

EFFECT OF LARGE CONVECTION HEATING 

In determining the heat transfer coefficients given in table I, the en- 
ergy equation and boundary conditions in reference 8 were taken as 


u, 3T = C^l 
ez 8z 2 

(8) 

: = 0 T = T 

P 

(9) 

■ 5 T ■ T sat 

(10) 


where a flat uniform vapor gap thickness 5 has been assumed (see fig. 2). 
The radial convection term in equation (8) has been assumed small and 
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thus neglected. The temperature T p is the supporting plate temperature 
and T ga j. is the saturation temperature of the liquid drop. The convective 
velocity is directed downward from the surface of the liquid drop. It 
was found by a solution of the momentum equations to be of the form 


“ = ---<Pl-V £ I 

3 A 2 M 


* 2 3 

5z - z 


) 


(ID 


where A is the cross sectional area of a drop of volume V. 

Equations (8) to (11) are incomplete in that the vapor gap thickness 6 
is not specified. However, an energy balance at the interface yield an ad- 
ditional equation. 

( 12 ) 

z=6 


-P v xw( 6) = -k 


3T 

3z 


In equation (12), the right hand side represents the heat conducted to the 
interface while the left hand side represents the equivalent latent heat re- 
moved from the interface. 

Equations (8) to (11) can be solved to yield the following expression 
for the temperature distribution and gradient 



(13) 


( 14 ) 


where a dimensionless gap thickness is defined as 
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with 6 * 
yields 


0 = 


2ttV( Pt 


P v )g 


L 3A Z a ii J 


6 4 = i'l 


r*\/C AT\ 

H 5 


*4 


*2 , 


(15) 


defined in table I. Substituting equation (14) into equation (12) 



Thus, the dimensionless gap thickness 0 is only a function of the param- 
eter (X/Cp AT). 

The heat transfer coefficient to the drop h is defined at the liquid va- 
por interface as follows: 


h < T p - TV = 9 = -*>/"(« 


(17) 


Substituting in the expression for 00 ( 5 ) from equation (11) into equation (17), 
introducing the parameter 0 from equation (15) and expressing the V and 
A in terms of V* and A* in table I give the expressions for the heat 
transfer coefficients shown in table I with X* in the dimensionless property 
group defined as 


X* = X 


,c p Af 



(18) 


where 0 can be found from a solution of equation (16) for various values of 
C AT/X as the parameter. 

r 

Exact Solution for X* 


Equation (16) was solved numerically to determine the value of 0 for 
various values of (C AT A)- A half-interval search technique in conjunc- 
tion with a Simpson rule integration was used. The value of 0 was then 
substituted into equation (18) to determine X*. Since the dimensionless va- 
porization time t* as well as the heat transfer coefficient are proportional 
to (X*) 1/4 , the fourth root of X*A is plotted in figure 3 instead of X*A- 
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Small C AT/X 
P 

For small values of C AT/x, the exponential in equation (16) can be 

tr 

approximated by 

e' x = 1 - x (19) 

Using this approximation in equation (16), solving for / 3 , and substituting 
the value for j 8 into equation (18) yields 


X small = X 


1 + 


C (T 
PIP 


T sat)l 


20 


( 20 ) 


This was the value used in the correlation of the data in figure 1. The 
fourth root of X* ma ^/X is also plotted in figure 3. 

Maximum X* 


The integral in the denominator of equation (16) cannot be evaluated in 
closed form. Consequently, the numerical and exponential approximations 
were necessary. However, another interesting solution for X* can be ob- 
tained by noting that 


■Kv 3 ~V A /2) 


-0T7/2 


e r-M ■, > g 0 < 7] < 1 

The expressions are equal at 77 equals 0 and 1. Therefore, 

r/3/2 


,C p AT 


w< 


X 


1 o - 0 t ?/2 


( 21 ) 


( 21 ) 


drj 


Solving equation (21) for (3 and substituting the expression for /3 into 
equation (18) yield 


X* <X| 


C AT' 
In / 1 + — 2 — 
2X 

C p AT 


= x* 

max 


( 22 ) 


2X 
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The fourth root of x* /\ is also plotted in figure 3. As seen in this 

max 

figure, X „ /X is always greater than the exact solution, 
max 

Surprisingly, however, the ratio of the fourth root of the exact to 
maximum solution is nearly a constant over a very large range of Cp AT A. 
This is seen in figure 3 by noting that the curves for X (exact) and * max 
are parallel. Consequently, the following slight modification of X max 
(eq. (22)) is made to give a good curve fit to the exact solution when 
C AT A is large: 

r 




f c at> 

0. 874 In fl + -P— 

2X 


C P AT 

2X 


C AT 

for — 2 > 2 


(23) 


As can be seen in figure 3, for large values of the parameter Cp ATA, 
the effect of convection heating becomes very important. The approximation 
Xgmaii should not be used in this case, rather, the numerical results or the 
curve fit given by equation (23) should be used. The Cp AT A will increase 
at higher system pressures, especially near the thermodynamic critical 
point. 

The nitrogen data in figure 1 has a Cp AT A value near 1. In this 
range the approximation for X* by X* ma ^ is still quite good. Conse- 
quently, the nitrogen data in figure 1 when corrected for the exact value of 
X* moved upward by only about 3 percent. This is an improvement in the 
theory; however, the correction is still small compared to the 30 percent 
deviation between the theoretical and experimental data. 


DISCUSSION 


Vapor breakthrough has been shown to be relatively small for the pres- 
ent set of data for ordinary fluids as shown in figure 1, although vapor break- 
through may still account for some of the downward shift in the experimental 
points for liquid nitrogen. A new theoretical curve which accounts for vapor 
breakthrough (shown dotted in fig. I) has been developed that seems to cor- 
relate the nitrogen data. 
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Convection heating effects, represented by the X* correction term, 
are relatively small for the temperatures and pressure (atmospheric) in- 
vestigated in this experiment. However, convection effects are likely to 
be more significant for higher temperatures and pressures near the criti- 
cal point. 

At present, the precise reasons for the downward shift in the nitrogen 
data in figure 1 has not been established. The most likely reason, at pres- 
ent, appears to be surface tension effects on the vapor dome. Water with 
a high value of surface tension tends to prevent the vapor from breaking 
through the domes, while nitrogen with a small value of surface tension al- 
lows the vapor domes to spew out the vapor. Thus, the flat disk model 
(upper curve in fig. 1) best approximates water and high surface tension 
fluid, while the breakthrough model (lower dotted curve in fig. 1) best ap- 
proximates nitrogen and lower surface tension fluids. 

In order to fully explain the anomalous results observed in figure 1, 
investigation of other possible causes for such behavior is warranted. 

Three such areas worthy of further investigation are surface tension ef- 
fects, mass transfer from the top of the drop and slip at the liquid vapor 
interface. The no slip boundary condition assumes the radial component 
of vapor velocity u (see fig. 2) is zero at the liquid vapor interface. This 
assumption works well for ordinary fluids. However, perhaps the rela- 
tively large difference in surface tension between ordinary fluids and the 
cryogen could alter the boundary condition. Wachters (ref. 3) shows that 
for the complete slip boundary condition the heat transfer coefficient in- 
creases by 40 percent. Consequently, a change in the boundary condition 
from slip to partial slip would improve the correlation, although justifi- 
cation for this partial slip boundary condition has not been physically estab- 
lished. 

As part of the present study, data was taken for large masses of water 
with surface tension additives. The data were correlated onto the plot using 
the value of the surface tension data without additives. If the correct value 
of a were known, the data point would shift to the right and upward. The 
strong downward trend in the water data point at a dimensionless volume 
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6000, however, indicates that perhaps surface tension effects should be 
investigated in more detail. 

Presently, the reason for the separation between the data for ordi- 
nary fluids and the cryogen remains unknown, in the sense that a variety 
of causes could produce the effect. 

CONCLUSIONS 


The principle conclusions of the present study are 

1. For large values of the parameter Cp AT/X the following expres- 
sion should be used for the modified latent heat of vaporization 


X cr= X 


, c at' 

0. 874 In ( 1 + J2 — 
2 \ 


2A 


C AT 

for -2 > 2 


This expression would also be used for conventional pool film boiling 
problems. 

2. For very large liquid nitrogen drops in film boiling, the following 
new expression has been found to correlate the vaporization time: 


r = 4.5 In i 


V* 

,14.2 


for V > 256 


3. Large differences in surface tension could have pronounced effects 
on the film boiling mechanism. 
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Figure 1. - Total vaporization time correlation. 
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Figure 2 . - Schematic model of the evaporation of a flat disk. 
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Figure 3. - Comparison of modified latent heats of vapor- 
ization. 


